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Abstract. We prove the convergence of the solutions of the parabolic wave equation to that of 
the Gaussian white-noise model widely used in the physical literature. The random medium is 
isotropic and is assumed to have integrable correlation coefficient in the propagation direction. We 
discuss the limits of vanishing inner scale and divergent outer scale of the turbulent medium. 



1. Introduction 



The small-scale refractive index variations, called the refractive turbulence, in the atmosphere is 
the result of small scale fluctuations of temperature, pressure and humidity caused by the turbu- 
lence of air velocities. For optical propagation in the atmosphere the influence of the temperature 
variations on the refractive index field is dominant whereas in the microwave range, the effect of 
the humidity variations is more important. The refractive turbulence results in the phenomena of 
beam wander, beam broadening and intensity fluctuation (scintillation). It is important to note 
that these effects depend on the length scales of the waves as well as the refractive turbulence |19j . 

The refractive turbulence is modeled on the basis of Kolmogorov theory of turbulence which 
introduces the notion of the inertial range bounded by the outer scale Lq (of the order of 100m — 1km) 
and the inner scale Iq (of the order of 1 — 10mm). Other features of the refractive turbulence in the 
open clear atmosphere include |22| : (i) small changes (typical value of 3 x 10~ 4 at sea level) in the 
refractive index related to small variations in temperature (on the order of 0.1 — 1°C), (ii) small 
scattering angle which is of the order and has the typical value 3 x 10 _4 rad for A = 0.6mn 
and £q = 2mm. Perturbation methods for solving the Maxwell equations are adequate provided 
that the propagation distance is less than, say, 100m, a severe limitation on their applicability to 
imaging or communication problems. Our motivation is mainly from laser or microwave beams but 
our consideration and results apply equally well to ultrasound waves in atmospheric turbulence. 
The results are also relevant in the context of ultrasound waves penetrating through complicated 
multiscale fluctuating (interface) zones in for instance human tissue. 

Under the condition A = 0(£q) (including the millimeter and the sub-millimeter range) the 
depolarization term in the Helmholtz equation for the electric field is negligible [221 and one can 
use the (scalar) Helmholtz equation 



with appropriate boundary conditions where k is the wavenumber, n is the mean refractive index 
field and h is the normalized fluctuation of the refractive index. 

1.1. The rescaled parabolic approximation. The well-known parabolic approximation to equa- 
tion Q is applicable in a regime where the variations of the index of refraction are small on the 
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scale of the wavelength so that backscattering is negligible P2J. This is almost always valid for 
laser beam in the atmosphere. 

In this paper we study the initial value problem for the parabolic wave equation 

(2) Vi*(z, x) + 2ik 9 ^^ X) = -2k 2 n(z, x)fr(z, x), *(0, x) = F (-) G L 2 (R 2 ) 



dz va, 
where z is the longitudinal coordinate in the direction of the propagation, x = {x\ } X2) is the 
transverse coordinates, Vj_ is the transverse gradient and \& is related to the scalar wave field E by 
E = ty(z,x) exp {iknz). The initial condition has a typical width a which is the aperture. Below 
we will drop the perp in denoting the derivatives in the transverse directions. 

The difficulty in solving equation © lies in the random multiscale nature of h(z,x). First we 
non-dimensionalize eq. (J2J) as follows. Let L z be the propagation distance in the longitudinal 
direction. Let Ao be the characteristic wavelength. The corresponding central wavenumber is 
&o = 2ti/\q. The Fresnel length Lf is then given by 

Lf = \/L z /ko. 

We introduce dimensionless wave number and coordinates 



k = k/ko, x = x/Lj, z = z/L 
and rewrite the equation in the form 

(3) 2ik— + + 2k 2 k L z n(zL z ,xL f )^ = 0, *(0,x)=. 

after dropping the tilde in the coordinate variables where 



T \ 2 



7 = 

is assumed to be 0(1), thus the source is supported on the scale determined by the Fresnel length. 

1.2. Model spectra. A widely used model for the structure function of the refractive index field 
of the atmosphere is based on the Kolmogorov theory of turbulence and has the following modified 
Von K' arm 'an spectral density 

(4) $ n (k) = 0.033C2(|k| 2 + ^ 2 )- 11 /6 e xp(-|k| 2 /O 

where k = (£, k), with £ G R, k G M? the Fourier variables conjugate to the longitudinal and 
transversal coordinates, respectively. Here Kq = 2ir j 'Lo,K m = 5.92/£o- This spectrum has the 
correct behavior only in the inertial subrange, i.e. 

(5) $ n (k) ~ |k|- 11/3 , |k| G (2ttLo 1 ,2^ 1 ). 

Outside of this range, particularly for |k| 2ttLq 1 there is no physical basis for their behavior; they 
are just mathematically convenient expressions of the cutoffs. In particular, if the wave statistics 
strongly depend on £q or Lq, then the problem probably requires more accurate information on the 
refractive index field outside of the inertial range [Sj, ^2]) Note that the ratio Lq/£q grows 

like Re 3 / 4 as the Reynolds number Re tends to infinity. 

There are several variants of @ arising from modeling more detailed features of the refractive 
index field. One of them is the Hill spectrum [2j, to account for the "bump" at high wave 
numbers which is known to occur near the inner scale 

(6£ n (k) = 0.033C 2 [l + 1.802|k|/K m - 0.254(|k|/^ m ) 7 / 6 l (|k| 2 + K 2 ) exp (-\k\ 2 /Kl) 



where K m = 3.3/ £q. The coefficient C 2 is itself a random variable that depends on time as well as 
the altitude. Note that in atmospheric turbulence the inner and outer scales and the exponent in 
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the power law may also have to be modeled as stochastic processes The temporal dependence 
is irrelevant for optical propagation; the altitude dependence has a rather permanent, non-universal 
structure with length scales much greater than the outer scale Lq |19j . 
We will consider a class of spectra satisfying the upper bound 

(7) *{hS) < K(L 2 + \k\ 2 y H - 3 / 2 (i + f \k\ 2 y\k = (£,k) eR 3 ,H e (0,1) 

with some constant K < oo as the ratio Lq/Iq — > oo in the high Reynolds number limit. The details 
of the spectrum are not pertinent to our results, only the exponent H is. In particular, H = 1/3 
for the modified Von K'arm'an spectrum (j3J). 

1.3. White noise scaling. Let us introduce the non-dimensional parameters that are pertinent 
to our scaling: 

L f L f L f 

£= \T ' V= T> P= f- 

In terms of the parameters and the power-law spectrum in ((7J) we rewrite (j5J) as 

f)\Il e h 2 7 , 

( 8 ) 2ik^- + ^ £ + — <7tfVK,x)tf e = 0, f £ (0,x)=F (7 1/2 x)e£ 2 (» 2 ). 
with 

Lf 

(9) a H = 

where p is the standard deviation of the refractive index field corresponding to <&(H, k). The 
spectrum for the (normalized) process V is given by 

(10) $ r/iP (k) < K( V 2 + \k\ 2 )- H - 3 / 2 (l + p- 2 |k| 2 )~ 2 ,k= (e,k) £R\H£ (0,1) 

which is rescaled version of l|7|). For high Reynolds number one has Lq/Iq = p/?) > 1 which is 
always the case in our study. 

In the beam approximation one has e<l. The beam approximation is well within the range of 
validity of the parabolic approximation. The white- noise scaling then corresponds to an = Oil). 
We set it to unity by absorbing the constant into V. This implies relatively weak fluctuations of 
the index field, i.e. 

C n ~ L 3 / 2 ~ H L-V 2 «1, as L z -> oo 

in view of the fact that H G (0, 1) and e < 1. 

In the present paper we first study the white-noise scaling with p < oo and rj > fixed as e — > 0. 
We then discuss the resulting white- noise model with p — > oo and rj — > 0. For the proof, we adopt 
the approach of ^01 where the turbulent transport of passive scalars is studied. In the white 
noise limit is studied via the so called Wigner distribution. 

2. Formulation and main results 
2.1. Martingale formulation. We consider the weak formulation of the equation: 

(11) ik[m,0)-{*o,0)] = -J Z ±(*l,A9)ds-^J Z (*%V(^,-)-9)ds 

for any test function 9 € C^°(1R 2 ), the space of smooth functions with compact support. The 
tightness result (Section 4.1) implies that for L 2 initial data the limiting measure P is supported in 
the Skorohod space D([0, zq]; L^QR 2 )). Here and below L^M 2 ) denotes the standard L 2 -function 
space with the weak topology. 
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For tightness as well as identification of the limit, the following infinitesimal operator A £ will 
play an important role. Let Vf = V{z/e 2 , •), T £ the cr-algebras generated by {Vf , s < z} and W z 
the corresponding conditional expectation w.r.t. T £ . Let A4 £ be the space of measurable functions 
adapted to {^J, Vz} such that sup^ <ZQ E|/(z)| < oo. We say /(•) £ £>(.4 e ), the domain of A £ , and 
A £ f = gif f, g eM £ and for = 5~ l [WJ{z + 8)- f(z)] we have 



su P E|/ d (z) 
KmE\f s (z)-g(z) 

o— »U 



< oo 

= 0, Vz. 



Consider the special class of admissible functions f(z) = cp((^ £ z ,9)), f'(z) = (ft'((^ £ z ,6)),V(j) £ 
C°°(K), then we have the following expression from (jllj) and the chain rule 



(12) A £ f(z) 
A main property of A £ is that 



-L(* £ Z ,A6) + -{* £ Z ,V £ Z 

2k £ 



(13) 

Also, 

(14) 



A £ f(s)ds is a .Ff -martingale, V/ £ V(A £ 



nm - f(s) 



E £ s A £ f(T)dr Vs<z a.s. 



(see HZ!). We denote by A the infinitesimal operator corresponding to the unsealed process V z {-) = 
V(z,-). 
Define 

J J cos((x-y).p)cos(aO*(^)(e,P)dadedp 
= vr y cos ((x - y) • p)^(^, p ) (0, p) dp 

(16) r«(x) = r«(x,x) 

where we have written the wavevector k £ R 3 as k = (£, p) with p £ M 2 . 

Now we formulate the solutions for the Gaussian Markovian model as the solutions to the cor- 
responding martingale problem: Find a measure P (of * z ) on the subspace of D([0, oo); L 2 (R 2 )) 
whose elements have the initial condition i ? o(7 1 ^ 2 x) such that 



2k 



is a martingale w. r. t. the filtration of a cylindrical Wiener process, for each f £ C c 



where 
(17) 



ic*.e = J ^(x)^( y )r( 1 )(x, y )0( y )d y . 



The Gaussian Markovian model has been extensively studied for beam wander, broadening and 
scintillation effects in the literature (see, e.g. ^1)- ^ can a ^ so been written as the Ito's 

l 



equation 



d9 z = ( -4A -Pr^ J ^ z dz + ik{K^ z ) l/2 dW z 
2k 

i 

2k 



-4 a - Pr^ ) V z dz + iJfetfadWi 



(18) = -^A^ z dz + ik^ z odW z , * (x) = F (7 1/2 x) 

where o stands for the Stratonovich integral, and W 2 (x) and W z (x) are the Brownian fields with 
the spatial covariance <5(x — y) and rW(x, y), respectively. 

The existence and uniqueness for the Schrodinger-Ito eq. (|18|) with p < oo and rj > has been 
studied in jS] by using the Wiener chaos expansion. Note that the following limit exists 

(19) T(x,y) = lim tt / cos((x-y) • p)$ (7?p) (0,p) dp 

= vrycos((x-y) • p)$ (r?)OO )(0,p) dp, 7] > 

By the well-posedness result of j5] and a standard weak—* compactness argument one can prove 
the existence of weak solution in -D([0, oo); L^(M 2 )) for p = oo, H € (0, 1). 

Next we consider the limiting case r] = 0. This would induce uncontrollable large scale fluctuation 
the Gaussian, Markovian model which should be factored out first. Thus we consider the solution 
of the form 

#(2,x) = y'(z,x)exp(ik / W s (0)ds) 

Jo 

and the resulting equation 



(20) JM Z = ^M> z dz + ik^ z odW' z , * (x) = F (7 1/2 x) 

2k 

where W' z is given by 

(21) W' z {x) = W z (x)-W z (0) 
with the covariance function 

F(x,y) = tt y (e-P - l)(e-*-P - l)<J> ( o !0o) (0,p)dp. 

Note that the above integral is convergent only if 

H < 1/2; 

in particular, the limit exists for the modified Von K'arm'an spectrum H = 1/3. Since H < 1/2, 
the limiting model is only Holder continuous in the transverse coordinates. 

Again by the well-posedness result of [Sj and a standard weak—* compactness argument one can 
prove the existence of weak solution in D([0, oo); L^IR 2 )) for p = oo, 77 = 0,iJ G (0, 1/2). 

2.2. Uniqueness. Because of the non-smoothness (when p = 00) and the non-homogeneity (when 
rj = 0) of the white-noise potential in the transverse coordinates the uniqueness argument of [HI 
does not apply here. 

Taking the function /(r) = r n in the martingale formulation, we arrive after some algebra at the 
following equation 

(22) ^j— = dF^+C 2 F^ 
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for the n— point correlation function 

Fj n )( Xl , . . . ,x n ) = E [tt,(xi) • • • * 2 (x n )] 

where 

n 

(23) Cl = 4^A Xj 

n 

(24) c 2 = -P r(x,> Xfc ), 

i,fc=l 

n 

(25) or C 2 = -k 2 ^ r '(xj,Xfc) 
(26) 

We will now establish the uniqueness for eq. (|22|) with the initial data 

F {n) (x l5 . . . , x ra ) = E [* (xi) • • • *o(x n )] , *o 6 ^ 2 (M 2 ). 

In the former case (|24|) C 2 is a bounded, Holder continuous function and we rewrite eq. (|22l) in 
the mild formulation 



= exp (zCi)F (n) + / 2 exp [(* - s)Ci]C 2 Fj n ) 



whose local existence and uniqueness can be easily established by straightforward application of 
the contraction mapping principle. By linearity, local well-posedness can be extended to global 
well-posedness. 

In the latter case (|25|l C 2 is unbounded, Holder continuous function with sub-Lipschitz growth. 
We first note that C 2 is non-positive everywhere since 

n „ 

J2 r'^Xfc) = ttJ 5>**» - l)^(e—P - l)d> (0>oo) (p)dp > 0. 

j,fc=i j fc 

Hence both C\ and C 2 are generators of one-parameter contraction semigroups on L 2 (M. 2n ), thus by 
the product formula (Theorem 3.30, |JJ) we have 



lim 

m— >oo 



Z Z 1 m 

exp (— Ci) exp (— C 2 ) i 7 = exp [z{d + C 2 )]F 
m m J 



for all -F S L 2 (R 2n ), which then gives rise to a unique semigroup on L 2 (M 2n ). 

2.3. Main assumptions and theorem. Let V z be a z-stationary, x-homogeneous square-integrable 
process whose spectral density satisfies the upper bound (|1U(1 . 

Let and .F+ be the sigma-algebras generated by {V s : Vs < z} and {Vs : Vs > z}, respectively. 
Define the correlation coefficient 

(26) p{t) = sup sup E [hg\ . 

E[hl=0,E[h 2 ] = l z +* 

E[g]=0,E[g2]=l 

Assumption 1. T/ie correlation coefficient p(t) is integrable 

When V z is a Gaussian process, the correlation coefficient /?(i) equals the linear correlation 
coefficient r(t) which has the following useful expression 

(27) r(t) = sup R(t - ti - T 2 ,k)g 1 (T 1 ,k)g 2 (T 2 ,'k)dkdT 1 dT 2 

91 ,92 J 



where 



R(t,k) = J e^*(^,k)de 



and the supremum is taken over all gi,g2 £ L 2 (M. d+1 ) which are supported on (—00, 0] x R d and 
satisfy the constraint 

(28) J R(t-t',k)g 1 (t, 1 k)gi(t',k)dtdt'dk = J R(t - t', k)g 2 (t, k)g 2 (t', k)dtdt'dk = 1. 

There are various criteria for the decay rate (e.g., exponential decay) of the linear correlation 
coefficients, see [IB] . 

Lemma 1. Assumption 1 implies that the random field 

/•oo 

V 2 (x) = / E 2 [Vt(x)]dt 
J 

is ^-homogeneous and has a finite second moment which satisfies the upper bound: 

roo poo 

E[V 2 Z ] < / / \E[E z [V t ]E z [V s ]]\dsdt 

J z J z 



( 2 

p(t)dt ' 



Proof. Consider in the definition of the correlation coefficient 

h l= E Z (V S ) 6L 2 (fl,P,f,) 

h 2 = v t a 2 (fi,p,f+). 

We then have 

|E [E,[V s (x)]E 2 [V t (x)]]| = |E [E 2 [V s (x)]V t (x)]| < p(t - z)E 1 / 2 [E 2 Z [V S ]] E 1 / 2 [V 2 ] . 

Hence by setting s = t first and the Cauchy-Schwartz inequality we have 

E[E 2 Z [V S ]] < p 2 (s-z)E[V 2 z ] 
|E[E 2 [V s (x)]E 2 [V t (x)]]| < p{t-z)p(s-z)E[V 2 z l s,t>z. 

Therefore 

f'OO f'OO 

E[V 2 ] < \E[E z [V t ]E z [V s ]]\dsdt 

J z J z 

(21) < E[V 2 ]Q°° p(t)dt\ 

which, together with the integrability of p(t), implies a finite second order moment of V 2 
Corollary 1. For each L, zq < 00 and p < 00,77 > there exists a constant C such that 



sup E 

Z<Z Q 

x|<t 



AV A2 



< C 



□ 



for all H £ (0,1), A > 1. 



Proof. Analogous to ()21j) we have 

E[AVJ 2 < E[AV Z } 2 

A straightforward spectral calculation shows that 

E[AV 2 (x)] 2 = 0(p 4 ~ 2H ) 



p(t)dt 



Vx, z. 



It is easy to see that 
and that 



AV Z 



r (i) ( 



E 



-V z 



V,(x)V z (y) 



where is given by (|15|). 

Next we assume the following quasi-Gaussian property: 



Assumption 2. 

(21) 

(22) 

(23) sup E 

\y\<L 



sup E[V/(y)] 4 < C\ sup E 2 [V/] 2 (y 



|y|<£ 

sup E 

|y|<£ 



Vf 



|y|<i| 

(y) < C 2 sup E 2 



Vf 



(y) < C 3 supE[V; £ ] 2 (y) 

\|y|<i 



+ sup E 2 

\\y\<L 



y £ y £ 

v z ' z 



(y) 









I sup E 2 




2 (y)j 


\lyl<^ 












I I sup E 


v( 


2 (y)J 


/ \\y\<L 







for all L < oo where the constants C\,Ci and C3 are independent of e, rj, p, 7. 
Assumption 3. For any fixed r\ > and every 6 G C^°(M 2 ) 



(23) 



SUp ||0V(-o,-)l|2 



o I - 



z<zo 



Ve < 1 < p 



with a random constant of finite moments independent of p ande. 

When V is Gaussian, V is also Gaussian and condition Q23|) is always satisfied 



(24) 



sup ||^V(_,.)|| 2 < Clog 



z<zo 



20 



□ 



where the random constant C has a Gaussian-like tail by a simple application of Borell's inequality 
&■ 

Theorem 1. Let V satisfy Assumptions 1, 2 and 3. Let n > and p < 00 be fixed as e — ► 0. Then 
the weak solution ^ £ of ill]) converges in the space of D([0, 00); L 2 U (M?)) to that of the Gaussian 

white-noise model with the covariance functions F^ and Iq . 

Note that in the limiting model with p = 00 the white-noise velocity field has transverse regularity 
of Holder exponent H + 1/2. 

The convergence of the white-noise limit has been established in and [3] for a refractive index 
field that is a function of z only. 



3. Proof of Theorem 1 
3.1. Tightness. In the sequel we will adopt the following notation 

m=f{(9i,o)), f(z)^fm,6)), f"(z)=rm,e)), v/gc c 

Namely, the prime stands for the differentiation w.r.t. the original argument (not z) of /, /' etc. 

A family of processes {^ e ,0 < e < 1} C D([0, oo); L^(]R 2 )) is tight if and only if the family 
of processes {(* £ , 6), < e < 1} C D([0, oo); L 2 W (R 2 )) is tight for all 9 G C C °°(R 2 ). We use the 
tightness criterion of 18 (Chap. 3, Theorem 4), namely, we will prove: Firstly, 



lim limsupP{sup | (¥|,0) | > N} = 0, Vz < oo. 

N~*oo £ ^Q Z<ZQ 



(25) 



Secondly, for each / G C°°(R) there is a sequence f £ (z) G T>(A £ ) such that for each zq < oo 
{^l e / e (2;), 0<e<l,0<z< zq} is uniformly integrable and 



(26) 



lim P{ sup \f £ (z) - f({^ £ ,9))\ >5} = 0, V5>0. 



z<zo 



Then it follows that the laws of {{$> £ , 9) , < e < 1} are tight in the space of D([0, oo); L 2 (M 2 )) 
Condition ()25|) is satisfied because the L 2 -norm is preserved. Let 



be the 1-st perturbation of /(#). Let 



E| /'(*) (%,V £ S 9) ds 



E £ z V £ ds. 



We obtain 
(26) 

Proposition 1. 



ft(z)=ikef(z)(%,V : 



lim sup E|/f (z)\ = 0, lim sup |/i(^)| = in probability 



Z<ZQ 



z<z 



Proof. We have 

(27) 

and 

(28) 



E[\fl(z)\] <e||/lUI|*o|| 2 E||0Vf||2 



sup 1/^(^)1 < ^H/'Hooll^olb sup \\9V £ Z 



z<z 



z<z 



The right side of l|27jl is 0(e) by Lemmanwhile that of (|28[) is o(l) in probability by Assumption 3. 

□ 

Set f £ (z) = f(z) + f £ (z). A straightforward calculation yields 



A £ ff 



-ef(z) 
-ef'(z) 



k 2 



(Vl,A9) + -(*l,V 



m,V £ z 9 



I(^,A(^)) + ^(*:.V:V^ 



ik 



f\z)^%m 



and, hence 
(27) A £ f(z) 



Yk f\z) (9l,AB) - k 2 f'(z) (n,V £ z V £ z e) ~ ~k 2 f"(z) (*:\ V:«) (*:. V: 
S \f{z) fa A(VI9)) + f"{z) (91, A9) fa VIO \ 



2 L 



A\{z) + Al{z) + A%{z)+Al{z) 



where A^(z) and A^(z) are the 0(1) statistical coupling terms. 

For the tightness criterion stated in the beginnings of the section, it remains to show 

Proposition 2. {A £ f £ } are uniformly integrable and 

lim supE|A£(z)| = 

e^o z<Zo 



Proof. We show that {Af }, i = 1, 2, 3, 4 are uniformly integrable. To see this, we have the following 
estimates. 

\A\{z)\ < i||/'|U|* o ||2||A0|| 2 
Ik 

\A £ 2 {z)\ < P||/'|UI*o|| 2 ||VfV^|| 2 

1^1(2)1 < fc 2 ll/"llooll*olllll^|| 2 ||v^e|| 2 . 

For fixed rj, the second moments of the right hand side of the above expressions are uniformly 
bounded as e — > and hence Af (z), Af,(z), A|(z) are uniformly integrable. 

; Viiooii^oii^A^ibiiv^ib + u/'iura2iiA(v^ 



< - 
~ 2 



By Lemma ^ and Corollary ^A| is uniformly integrable. Finally, it is clear that 

lim supE|AJ(z)| = 0. 



z<z 



□ 



3.2. Identification of the limit. Once the tightness is established we can use another result in 
|18j (Chapter 3, Theorem 2) to identify the limit. The setting there is finite-dimensional but the 
argument is entirely applicable to the infinite-dimensional setting here (cf. jllj). 

Let A be a diffusion or jump diffusion operator such that there is a unique solution uj z in the 
space D([0,oo);Ll (R 2 )) such that 

(22) f(u z )- f Af(oJ s )ds 

Jo 

is a martingale. We shall show that for each / E C°°(R) there exists f e G T>(A £ ) such that 
(23) 

(24) 
(25) 

(26) 



sup E\f £ (z)-f((9 £ z ,6))\ < oo 

Z<Z0,£ 

0, \/z < Zq 



limE|/ 6 (z) -/(<¥£, 
sup E\A £ f(z)-Af((9l,e))\ < oo 

Z<ZQ,£ 



\imE\A e r(z)-Af((9l,e))\ 
e— >0 



0, \/z < Zq. 



Then it follows that any tight processes {9 £ ,9) converges in law to the unique process generated 
by A. As before we adopt the notation f(z) = f((9 £ z ,0)). 
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For this purpose, we introduce the next perturbations /|,/|- Let 
(27) A^(cf>) = y , y^(x)0(x)r( 1 )(x,y) ( A(y)^(y) dxdy 



(28) 

where 
(29) 

It is easy to see that 
(30) 



A 



= J rW(x,x)0(x)0(x) dx 
r( 1 )(x,y) ^ E[v|(x)V|( y ) 
4 1) (0)=Efev^)(<A,V^ 



Define 



ds 



fl(z) ee k 2 f'(z 



Ef 



ds. 



Let 

and 
(31) 
(32) 

we then have 
(33) 

(34) 

Proposition 3. 



r( 2 )(x,y)^Efe(x)V z £ (y)' 



A ?\<P) = j y^x)0(x)r( 2 )(x,y)0(y)e(y) dxdy 
4 2) (0) = J r( 2 )(x,x)<Xx)fl(x) dx, 



2 1,2 



/!(*) 



e'k 
e 2 k 2 



lim sup E|/f (z)| = 0, lim sup E|/f (z)| = 0. 



z<z 



z<z 



Proof. We have the bounds 



supE|/f(z)| < sup^^HnU |^o||^E||V^||^ + E[4 2) (^)] 

z<zq z<zq L 



z<z 

both of them tend to zero. 



supE|/|(z)| < supe 2 k 



2 7,2 II t i 



Z<ZQ 



|^ ||2E||V|V^|| 2 + E[4 2) (^)] 



We have 



A £ fl(z) 



k 2 f"(z) -(n,VP}(n,VP) + AX>(n) + 



li 



with 



m(z) 



ie k 



27f'"(z) 



+u?kf''(z)(* e z ,v e g e 



(29) 



-ie 2 ~kf"(z) 



(2) 

where G e denotes the operator 

Gf<t> = f e(x)r( 2 )(x,y)^(y)0(y) ( iy. 



Similarly 



ie 2 kf'(z] 
+ie 2 kf"(z) 
-ie 2 kf'(z) 



i(^,A(Vf^)) + ^(*:.V:V:V:« 



4\n 



^(n,A(v^9)) + ^(%,v^t 



where 



Proposition 4. 



if(x)^r( 2 )(x,x). 



lim sup E\R £ 2 (z)\ = 0, lim sup E\R £ 3 (z)\ = 0. 

£->0 2 < zo £^0 z<zo 



The argument is entirely analogous to that for Proposition |21 The most severe factors involve 
A(VJ0) and A(V £ V £ 0), both of which have uniformly bounded second moments by Assumption 2 
and Corollary ^ Therefore the corresponding terms are 0(e 2 ). 

Consider the test function f £ (z) = f(z) + f £ (z) — f 2 (z) — fli z )- We have 



(27) A £ f £ {z) 



= ±f(z) (%, A9) - k 2 f"(z)4\%) - Pf'A^m) - Bftz) - R £ 3 (z) + A\{z) 
2k 

Set 

(27) R £ (z) = R\ (z) - R £ 2 {z) - Rl (z) , with R\{z) = A\{z). 

It follows from Propositions 3 and 5 that 

lim supE|i? e (z)[ = 0. 



u2 



£^0 



Z<ZQ 
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Recall that 



Mm 



A £ f £ {s)ds 



f(z) + f!(z)-fl(z)-fl(z) 



-f-f'iz) (91, A9) ds 
Ik 



+ 



f"{s)A ( }\%) + f\s)A\ L \9 



ds 



R e (s)ds 



is a martingale. Now that (|23|) - (|26|) are satisfied we can identify the limiting martingale to be 



(24) 



M z {9) = f(z) 



f'(s) 



k 2 f"(s)A^'(9 s )}ds. 



i(i) 



Since {9 £ z ,9) is uniformly bounded 

K*!,0)l<ll*ol|2||0|| 2 

we have the convergence of the second moment 

limE{<**,0) 2 } =e{(*„i 

Use f(r) = r and r 2 in (|3.2|) 



'}■ 



M^(9) = (9 Z ,9)- 
is a martingale with the quadratic variation 



where 



Therefore, 



ds 



A { ^(9 s )ds 



AfW(0),AfW(0) 

1C* 3 9 = j # s (x)r«(x,y)^(y)0(y)dy. 

=ik I y/K^jMs 
Jo 



?, 1C^ 3 9) ds 



where W s is a real- valued, cylindrical Wiener process (i.e. dW z (x) is a space-time white noise field) 
and y / /C* s is the square-root of the positive-definite operator given in (|17|). 
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